Free flux flow resistivity in strongly overdoped high-T c cuprate; purely viscous motion 
of the vortices in semiclassical rf-wave superconductor 
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We report the free flux flow (FFF) resistivity associated with a purely viscous motion of the vor- 
tices in moderately clean d-wave superconductor Bi:2201 in the strongly overdoped regime (T c =16 K) 
for a wide range of the magnetic field in the vortex state. The FFF resistivity is obtained by mea- 
suring the microwave surface impedance at different microwave frequencies. It is found that the 
FFF resistivity is remarkably different from that of conventional s-wave superconductors. At low 
fields (H < 0.2H C 2) the FFF resistivity increases linearly with H with a coefficient which is far 
larger than that found in conventional s-wave superconductors. At higher fields, the FFF resistivity 
increases in proportion to y/TI up to H C 2- Based on these results, the energy dissipation mechanism 
associated with the viscous vortex motion in "semiclassical" d-wave superconductors with gap nodes 
is discussed. Two possible scenarios are put forth for these field dependence; the enhancement of 
the quasiparticle relaxation rate and the reduction of the number of the quasiparticles participating 
the energy dissipation in d-wave vortex state. 

PACS numbers: 74.25.Fy 74.25.Nf 74.60.Ec 
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I. INTRODUCTION 

When a vortex line in a type-II superconductor moves 
in the superfluid, the frictional force is determined by 
the damping viscosity, which in turn depends on the en- 
ergy dissipation processes of quasiparticles. The prob- 
lem of the energy dissipation associated with the vis- 
cous motion of the vortices has continued much atten- 
tion of researchers for years. To gain an understanding 
on the energy dissipation, the experimental determina- 
tion of the free flux flow (FFF) resistivity is particu- 
larly important. Hereafter the term FFF will refer to 
a purely viscous motion of the vortices, which is realized 
when the pinning effect on the vortices is negligible. The 
FFF resistivity is known to be one of the most funda- 
mental quantities in the superconducting state. In fully 
gapped s-wave superconductors, the flux flow state has 
been extensively studied and by now a rather good un- 
derstanding on the the energy dissipation processes has 
been achieved. Jl|, ||, ||, ||, ||, ||, |j In s-wave supercon- 
ductors, the quasiparticles trapped inside the vortex core 
play a key role in the dissipation processes. Moreover, it 
has been shown that there is a fundamental difference 
in the quasiparticle energy relaxation processes among 
dirty (£ > £), moderately clean (£ < I < £ • ^) and su- 
perclean (i > £ • s-wave superconductors, where £ is 



the coherence length, £ is the mean free path, sp is the 
Fermi energy, and A is the superconducting energy gap. 

A renewed interest in the problem concerning the 
quasiparticle dissipation owes to recent developments 
in the investigation of unconventional superconductors. 
The latter are characterized by superconducting gap 
structures which have nodes along certain crystal direc- 
tions. In the last two decades unconventional super- 
conductivity has been found in several heavy fermion, 
organic and oxide materials. From the viewpoint of 
the physical properties of the vortex state, perhaps 
the most relevant effect of the nodes are the existence 
of gapless quasiparticles extending outside the vortex 
core.||, |9[ [l0| In fact recent studies of heat capacityOl, 
thermal conductivity ]T^| , and NMR relaxation rate|13| 
provide a strong evidence that these quantities are gov- 
erned by delocalized quasiparticles. However, despite 
these extensive studies of the vortex state of unconven- 
tional superconductors, the microscopic mechanisms of 
the energy dissipation associated with the viscous vor- 
tex motion is still far from being completely understood, 
exposing explicitly our incomplete knowledge of vortex 
dynamics in type-II superconductors. Thus it is particu- 
larly important to clarify whether the arguments of the 
energy dissipation are sensitive to the symmetry of the 
pairing state, jlq 
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Recently, the flux flow resistivities in /-wave supercon- 
ductor UPt3 and ci-wave high-T c cuprates, both with line 
nodes, were demonstrated to be quite unusual. However, 
these materials may not be suitable for the study of the 
typical behavior of the flux flow resistivity in unconven- 
tional superconductors. The T vs H phase diagram of 
UPt3, which still is controversial, is considered to con- 
sist of various phases with different superconducting gap 
functions, which complicates considerably the interpre- 
tation of the FFF resistivity. jlj], [ll| The flux flow resis- 
tivities of YBa 2 Cu 3 07_5 and P^Si^CaC^Og+s in the 
underdoped and optimally doped regimes have been mea- 
sured by several groups but here again there are several 
difficulties in interpreting them. 
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For instance, the measurements could not cover a wide 
field range in the vortex state due to extremely large up- 
per critical field H C 2- Moreover, very recent STM mea- 
surements have demonstrated that the vortex core struc- 
ture of these high-T c cuprates is very different from that 

tected in the semiclassical d-wave superconductor) 24 
possibly due to the extremely short coherence lengths 
and the strong antiferromagnetic fluctuation effect within 
the core. 

The situation therefore calls for the need for a text- 
book example of the FFF resistivity of unconventional 
superconductors with nodes, in which the semiclassical 
description of the vortex core discussed in the litera- 
ture, e.g. Refs.g, g, ffij |§ applies. Especially the FFF 
resistivity in the "semiclassical" superconductors in the 
moderately clean regime is strongly desired, because al- 
most all unconventional superconductors fall within this 
regime. It should be noted that the determination of the 
FFF resistivity is not only important for understanding 
the electronic structure in the vortex state but is also rel- 
evant for analyzing the collective motion of the vortices, 
such as the flux creep phenomena. This is easily under- 
stood if one recalls that the motion of the vortices in the 
vortex liquid and solid phase in high-T c cuprates has been 
analyzed by assuming the Bardeen-Stephen relation for 
individual vortex, as discussed in §V. Then, if the FFF 
resistivity strongly deviates from the Bardeen-Stephen 
relation, the interpretation of the collective motion of 
the vortices should be modified. 

We stress here that high-T c cuprates in strongly over- 
doped regime are particularly suitable for the above pur- 
pose because of the following reasons, (i) Most impor- 
tantly, it appears that the semiclassical description of 
the electronic structure of the vortex core is adequate 
in strongly overdoped materials. || ||, [n| This is 
because many experiments have revealed that in the over- 
doped regime the electron correlation and antiferromag- 
netic fluctuation effects, which might change the vortex 
core structure dramatically as observed in STM measure- 
ments, are much weaker than those in optimally doped 
and underdoped materials. In fact most of the physical 
properties in the overdoped materials are well explained 
within the framework of the Fermi liquid theory, (ii) 
Low H C 2 enables us to measure the FFF resistivity for a 



wide range of the vortex state, (iii) The large coherence 
lengths and small anisotropy ratio reduce the supercon- 
ducting fluctuation effect which make the interpretation 
of flux flow resistivity complicated. In fact, as we discuss 
in §IV, the resistive transition of the overdoped materials 
in magnetic field is much sharper than that of optimally 
doped and underdoped materials, (iv) The flux flow Hall 
angle which complicates the analysis of the flux flow state 
is very small. Q 

The purpose of this work is to present and discuss our 
experimental results on the FFF resistivity pj of mod- 
erately clean d-wave superconductors. The experiments 
were carried out using strongly overdoped Bi:2201. This 
system is an excellent choice for studying the FFF resis- 
tivity. It has a comparatively simple crystal structure (no 
chain, single Cu02 layer) and hence the band structure is 
simple. H C 2 is within laboratory reach over a very broad 
range of temperatures. A major cause of difficulty in ob- 
taining the FFF resistivity in high-T c cuprates was the 
strong pinning effect. To overcome this difficulty, we have 
measured the microwave surface impedance at different 
frequencies. High frequency methods are suitable for this 
purpose because they probe vortex response at very low 
currents when the vortices undergo reversible oscillations 
and they are less sensitive to the flux creep. [2^, We 
show that the FFF resistivity of the "semiclassical" d- 
wave superconductor is very different from that of con- 
ventional s-wave superconductors. On the basis of the 
results, we discuss the dissipation mechanism associated 
with viscous motion of the vortices in unconventional su- 
perconductors. 



II. EXPERIMENT 

High quality single crystals of Bi:2201 
(Bi1.74Pb0.38Sr1.88Cu1.00Oy) in the overdoped regime 
with transition temperature T c = 16 K were 
grown by the floating zone method. |2j| The sam- 
ple size used for the microwave measurement was 
^0.8mmx0.7mmx0.04 mm. The upper inset of Fig. 1 
depicts the magnetization at the superconducting 
transition for the same sample used for the microwave 
measurements. The normal state resistivity in the 
a6-plane p n depends on T as p n cx T 13 with (3 ~ 2; the 
typical Fermi liquid behavior which can be seen in the 
overdoped high-T c cuprates. The resistive transition 
of the sample in the same batch with T c =18 K is 
also shown in Fig. 1. Both resistive transition in zero 
field and magnetization measurements show a sharp 
superconducting transition. 

The microwave surface impedance Z s = R s + iX s , 
where R s and X s are the surface resistance and sur- 
face reactance, respectively, was measured by the stan- 
dard cavity perturbation technique using cylindrical cav- 
ity resonators made by oxygen free Copper operated in 
TEon mode. The resonance frequencies of these cavi- 
ties were approximately 15 GHz, 30 GHz, and 60 GHz. 
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T (K) 

FIG. 1: The resistive transition in magnetic field of overdoped 
Bi:2201 in the same batch with T c —18 K. Inset (upper): The 
magnetization at 5 Oe of the same sample used for the mi- 
crowave measurements under the conditions of zero field cool- 
ing (ZFC) and field cooling (FC). Inset (lower); T-dependence 
of H c 2 determined by three different methods. The filled tri- 
angles denote H C 2 defined by the dc-resistive transition in the 
main panel, using criteria p = 0.5p„. The filled circles denote 
H C 2 defined by the magnetic field at which p\ becomes fre- 
quency independent. The open squares denote K C 2 defined 
by the field at which R s reaches to a normal state value. H C 2 
is estimated to be ~ 20 T below 5 K. 

The sample was placed in an antinodes of the oscillatory 
magnetic field H ac , such that H ac lies parallel to the c- 
axis of the sample. The external dc-magnetic field was 
applied perpendicular to the afr-plane. In this configu- 
ration, the two dimensional pancake vortices respond to 
an oscillatory driving current induced by H ac within the 
a6-planes. The cavities at 15 GHz and 30 GHz were op- 
erated at 1.7 K and sample temperatures were controlled 
by hot finger techniques using sapphire rod. The sample 
temperature in the cavity at 60 GHz was controlled by 
changing the temperature of the cavity. The Q- values of 
each cavity are 6.2xl0 4 for 15 GHz, 2.3xl0 4 for 30 GHz 
at 4.2 K , and 2xl0 4 at 4.2 K and 1.5xl0 4 at 20 K for 
60 GHz. According to the cavity perturbation theory, R s 
and X s can be obtained by 

and 

X s =G(-^j-^)+C = G(-^j+C, (2) 

where Q s and f s are the Q-factor and the resonance fre- 
quency of the cavity in the presence of the sample, and 




T (K) 

FIG. 2: T-dependence of the surface resistance R s (a) and 
surface reactance X s (b) at 15 GHz in magnetic field. Both 
the microwave magnetic field H ac and dc magnetic field B 
are applied parallel to the c-axis (-H" a c|| B|| c). In this con- 
figuration, the energy dissipation is caused by the oscillation 
of the two-dimensional pancake vortices. The measurements 
have been done under the field cooling condition. The abso- 
lute value of R s and X 3 were determined by the normal state 
dc resistivity. Inset: AA = A(0) — A(T) at low temperatures 
is plotted as a function of T 2 . 



Qo and /o are those without sample. G is a geometrical 
factor and C is a metallic shift constant. 

In Figs. 2(a) and (b), the T-dependence of R s and 
X s for Bi:2201 at 15 GHz are shown. The measure- 
ments in magnetic field have been performed in the field 
cooling condition. We first discuss R s and X s in zero 
field. In zero field, both R s and X s decrease rapidly 
with decreasing T below the transition. Let us quickly 
recall the behavior of Z s in the superconductors. In the 
normal state, the microwave response is dissipative and 
R s = X s = pqujS, where po is the vacuum per meability, 
lu/2it is the microwave frequency, and S n = \j2p n j 'pqlo 
is the normal state skin depth. In Bi:2201, t is estimated 
to be ^200 A, which is well shorter than 6 n at the onset 
in our frequency range. We therefore can determine the 
absolute value of R s and X s from the comparison with 
p n assuming R s = X s (Hagen- Rubens relation). In the 
Meissner phase, the microwave response is purely reactive 
and R s ~ and X s = puj\ a i,, where A a /, is the London 
penetration depth in the a6-planc. Using p n = 130/iil 
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cm for Bi:2201 at the onset, we obtained A a b=1500 A at 
T = 0. This value is slightly smaller than the penetration 
depth in YBa 2 Cu307_5 and I^S^CaCi^Os+a. In the 
inset of Fig. 2(b), AA = A(0) — A(T) at low temperatures 
is plotted as a function of T 2 . A A is proportional to T 2 . 
The relation AA oc T 2 has been observed in many high-T c 
cuprates and discussed in terms of the supcrfluid den sity 
in d-wave superconductors with the impurity state. (3CJ 



III. SURFACE IMPEDANCE IN THE VORTEX 
STATE 

We now focus on the surface impedance in the vortex 
state. Figure 3 shows the iJ-dependence of R a and X a 
of Bi:2201 at 15 GHz. In these measurements R s and 
X s are obtained by sweeping H . The hysteresis due to 
the effect of the trapped field in the crystal is very small. 
Moreover, both R s and X s obtained by sweeping H well 
coincide with those obtained under the field cooling con- 
ditions shown in Fig. 2. These results indicate that nei- 
ther inhomogeneous field distribution inside the crystal 
nor magnetostriction |42) caused by sweeping H seriously 
influences the analysis of Z s . 

In the vortex state, Z s is governed by the vortex dy- 
namics. We may roughly estimate R s in the limit of large 
and negligible rf field penetration as follows. In the flux 
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flow state when the pinning frequency uj p /2ti is negli- 
gible compared to the microwave frequency (lu p <C u>), 
two characteristic length scales, namely A a b and the flux 
flow skin depth 5f ~ y/2pf / p^oj, appear in accordance 
with the microwave field penetration. At low fields, X a b 
greatly exceeds 5f (X a b Sf). In this regime, R s and 
X s are given as R s ~ Pf/X a b and X s ~ pou>X a b- On the 
other hand, at high fields where Sf greatly exceeds X a b 
(Sf 3> X a b), the viscous loss becomes dominant and the 
response is similar to the normal state (i? s ~ X s ) ex- 
cept that S n is replaced by Sf. In the presence of pinning 
centers of the vortices, R s is reduced as discussed below. 

We here analyze the field dependence of Z s in accor- 
dance with the theory of Coffey and Clem. |Q The equa- 
tion of vortex motion for the vortex line velocity it, 



-qu + k p x — $0 J x z 



(3) 



where rj and k p are the viscous drag constant and pinning 
parameter, respectively, and z the unit vector parallel to 
B (we take «7"|| a?) . According to Coffey and Clem, the 
field dependence of Z s in the Meissner and vortex phases 
is expressed as 



Z s = i^oUlXab 



1 - ^/2)S V 2 /X 2 ab 



-1 1/2 



1 



(4) 



where S 2 = 8Al — iuj p /uj)~ l with w p /27r = K p /2irr] being 
the pinning frequency. Writing Z s in ter ms of the com- 
plex resistivity p = p\ + ip2 as Z s = \J iujpo(pi + ip^), 
we have 



pi = p uj 



^ib s 
1 + s 2 



pi 



1 l + sp 

1 + s 2 1 + p 2 ' 



(5) 
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FIG. 3: Field dependence of the surface resistance R a (a) and 
surface reactance X s (b) at 15 GHz measured by sweeping H . 



FIG. 4: The field dependence of A Pl (H) = Pl (H) - p x (0) 
obtained from R a and X 3 at three different microwave fre- 
quencies, pi increases with increasing microwave frequency. 
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and 



p 2 = Po^ 



1 



Pf 



1+p 2 



(6) 



»/ 



and p = LOp/u. In Eqs.(]5|) and 



where s 

the first terms in the right hand side are pi and pi at zero 
field, and second terms represent the field dependence. In 
what follows we discuss the microwave response focusing 
on pi obtained from R s and X s . Figure 4 shows the field 
dependence of Api(H) — pi (if) — pi(0) at three differ- 
ent microwave frequencies. The field dependence of Api 
is frequency dependent; pi increases with increasing fre- 
quency. Since p\ is reduced by the vortex pinning effect, 
as seen in Eq.(^), this result indicates that the pinning 
effect of the vortices is not negligible for the analysis of 
the flux flow resistivity in our microwave frequency range. 
Therefore, it is necessary to determine the pinning fre- 
quency for an accurate determination of the FFF resis- 
tivity. 

In Fig. 5, Api at T=3 K is plotted as a function of the 
microwave frequency. The solid lines show the results of 
the fitting by Api(H,uj) = pfuj 2 /(uj 2 + lo 2 ). It should 
be noted that since s <JC 1 except the vicinity of H C 2, 
the -ff-dependence of s little influences the present anal- 
ysis. Nevertheless, we restrict our analysis at H <10 T 
to avoid the influence of ii-dependence of s. The fit- 
ting parameters are oj p and p/. The ambiguity for de- 
termining uj p and pf is small. The ii-dependence of the 
pinning frequency obtained by the fitting is depicted in 
Fig. 6. At low field, oj p /2tt is approximately 22 GHz 
at T =3 K and 17 GHz at 5 K. These values are much 
larger than the pinning frequency in Bi2Sr2CaCu2 08+<s 
but much smaller than w p /2-k in YBa2Cu307 [[l9| |2^, p3[ . 



At low field, uip decreases gradually, while at >1.5 T lu p 
decreases approximately as lo v cx H , as shown in the 
inset of Fig. 6. 



IV. FREE FLUX FLOW RESISTIVITY OF 
BL2201 

Before discussing the FFF resistivity, it will prove use- 
ful to first comment on H C 2 of Bi:2201. It is well known 
that the resistive transitions of high-T c cuprates are sig- 
nificantly broadened in magnetic field due to the strong 
thermal fluctuation effect and the vortex dynamics. Al- 
though in overdoped Bi:2201 such a broadening effect 
is relatively small, it still becomes an obstacle in deter- 
mining H c2 10- ln the lower inset of Fig. 1, we plot 
H C 2 determined by three different methods. The filled 
triangles represent H C 2 defined by the dc-resistive tran- 
sition in Fig. 1, using a criteria p — \p n - The filled 
circles are H C 2 defined by the magnetic field at which 
pi becomes frequency independent. The open squares 
represent H C 2 defined by the field at which R s reaches 
to a normal state value. The values of H c2 obtained 
from the three different methods do not differ signifi- 
cantly A striking divergence in H C 2 as the temperature 
approached zero was reported in the overdoped Tl:2201 
in the transport measurements, J32]| while such a diver- 
gent behavior was not observed in the specific heat and 
Raman scattering measurements. ]33| The divergent be- 
havior of H C 2 was discussed in terms of several proposed 
models, such as the Josephson coupled small grains with 
T c higher than the bulk.) 34, p3 However in the present 



Bi:2201 such anomalies are not observed in H C 2 at least 




FIG. 5: Frequency dependence of Api(if) at T =3.0 K. 
(Filled triangles (1.0 T), open squares (0.8 T), filled squares 
(0.6 T), open circles (0.4 T), filled circles (0.2 T)). The 
solid lines are the results of the fitting by Api(B,u) — 
pfuj 2 /(uj 2 + ujp). For details, see the text. 



FIG. 6: The field dependence of the pinning frequency u} p /2n 
at T = 3 K and 5 K obtained by the fitting shown in Fig. 4. 
Inset: Same data plotted as a function of 1/H. uj p /2n de- 
cays in proportion to 1/H. The sloid lines show the relation 
uip/2-K oc 1/H. 
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above 2 K. At present we do not know the reason for this 
difference. From these measurements, H c2 is estimated 
to be approximately 20 T below 5 K. 

In Fig. 7(a), we plot Pf/p n as a function of H/H C 2 
at 3 K, assuming H C 2=19 T. If we assume H C 2=17 T at 
T =5 K, both pf almost exactly coincide with pf at 3 K, 
as shown in Figs. 7(a) and (b). The field dependence of p/ 
is convex. We found that the there are two characteristic 
regimes in the 7?-dependence of pj. In the low field region 
(H/H C 2 <0.2), pj increases linearly with H as 



Pf 



H 

u-n—p-n 

tlc2 



(7) 



with a ~ 2. A deviation from ii-linear dependence is 
clearly observed at higher field. In Fig. 7(b), Pf/p n is 
plotted as a function of ^H/H C 2. We found that pj 
increases as 



Pf 



H 

H C 2 



(8) 



at H/H C 2 >0.2. Since the linear extrapolation of Pf/p n 
in Fig 7(b) points to Pf / p n = 1 at H / H C 2=1, it is natural 
to expect that the relation of Eq. m) continues all the way 
up to H C 2- 



DISCUSSION 



A. Flux flow in s-wave superconductors 

In order to contrast the present results with the FFF 
resistivity of isotropic s-wave superconductors, we first 
briefly review the flux flow state in s-wave superconduc- 
tors. 

For isotropic s-wave pairing in the dirty regime, the 
Bardeen-Stephen model appears to be quite successful in 
describing the energy dissipation, jl], |[ The Bardeen- 
Stephen theory models the vortex core as a cylinder 
whose radius is the coherence length. It is assumed that 
the core is a normal metallic state inside of which the 
energy dissipation is dominated by the impurity scat- 
tering, similar to the ordinary resistive process. This 
is a good approximation for dirty superconductors with 
I < £. It follows from this model that the FFF resis- 
tivity in dirty s-wave superconductors is proportional to 
the normal state resistivity and is to the number of the 
vortices, 



Pf = p n H/H c2 . 



(9) 



The validity of this Bardeen-Stephen relation has been 
confirmed in most of dirty s-wave superconductors al- 
most throughout the whole Abrikosov phase; H c i < H < 

Hc2- § 

However, the description of vortex core as a normal 
metal is limited to dirty s-wave superconductors. In 
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FIG. 7: (a) The flux flow resistivity at T =3 K and 5 K as a 
function of H/H c2 . We assumed H c2 =19 T at 3 K and 17 T 
at 5 K. The flux flow resistivity is normalized by the norma l 
state value, (b) Same data plotted as a function of H/H c2 . 



moderately clean and superclean s-wave superconduc- 
tors with i > £, the quasiparticle response to an elec- 
tromagnetic field is radically different from that of nor- 
mal electrons, since the model of a normal metallic core 
breaks down, g The differ ence lies in the fact that 
the quasiparticles in the core are subject to Andreev re- 
flections by the pair potential and form the bound states 
of Caroli, de Gennes and Matricon j3(], ^7) before get- 
ting scattered by impurities. The largest energy differ- 
ence between the bound states is roughly estimated as 
hflo ~ A 2 /ef, where fio is the angular velocity. The 
electric conduction in the vortex state is governed by the 
scattering time between the Andreev bound states in the 
presence of impurities. Effects of these quasiparticles on 
the vortex dynamics have been considered in a number 
of papers. For moderately clean s-wave superconductors, 
the FFF resistivity has been calculated asfl, M 



1 



H 



Pf ~ Pr, 



hi 



(10) 



The logarithmic factor results from the shrinkage of the 
vortex core at low temperature and logarithmic energy 
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dependence of the impurity scattering rate of the Andreef 
bound state (Kramer-Pesch effect) [B8| . Thus, in spite of 
the fundamental difference of the character of the quasi- 
particles within the vortex core, the FFF resistivity in 
the moderately clean s-wave superconductors increases 
in proportion to H, which is similar to that in the dirty 
superconductors. In fact, the FFF resistivity of several 
moderately clean s-wave superconductors discovered re- 
cently were found to be proportional to H, though the 
logarithmic correction at very low temperature has never 
been reported so far. p9| 



B. Flux flow in d-wave superconductor 

We are now in position to discuss the FFF resistivity 
of semiclassical <i-wave superconductors. It is obvious 
from Figs. 7(a) and (b) that the field dependence of pf 
expressed as Eqs. (j^J) and ^) is markedly different from 
that of conventional s-wave superconductors expressed as 
Eqs. (||) and 

We first discuss the low field behavior of the FFF re- 
sistivity in Bi:2201. The linear dependence of p/ on the 
magnetic field means that the energy dissipation per vor- 
tex does not depend on the magnetic field or the inter- 
vortex spacing. We can interpret this fact naturally if 
the energy dissipation is assumed to occur mainly near 
each vortex even in the superconductors with gap nodes. 
In fact, this assumption is justified by a numerical result 
on the ac response of the d-wave vortex. || Compar- 
ing Eq.(0) with Eq.([io|), the coefficient of the H- linear 
term in d-wave superconductor is found to be nearly as 
twice as that in s-wave superconductors. This behav- 
ior is similar to UPt3 with line nodes, in which pf at 
low field is larger than that found in conventional s-wave 
superconductors. 1 14 It should be noted that a similar 
result was reported in very recent measurements of high 
purity borocarbide superconductor YM2B2C with very 
anisotropic superconducting gap, presumably anisotropic 
s-wave symmetry. H^, These results lead us to con- 
clude that a large initial slope is a common feature in 
the FFF resistivity of the superconductors with nodes. 
In what follows, we discuss possible origins for the en- 
hancement of the FFF resistivity at low fields on the 
basis of the theoretical results available at the present 
stage. 

According to Kopnin and Volovik, the vortex trans- 
port in semiclassical d-wave superconductors is gov- 
erned by dynamics of quasiparticles which form An- 
dreev bound states around a vortex, much like in s-wave 
superconductors. |4l| The excitation spectrum of those 
quasiparticles is given by 



E(L,6) = -Q{9)L 



(11) 



in terms of the angle 9 in the momentum space and L 
the angular momentum. In this expression, 0,(6) denotes 
the angular velocity, which depends on the direction 9. 
Roughly speaking, fl(9) is proportional to the square of 



the energy gap, A(9) (oc cos(2#) for d x 2_ y 2 states). This 
branch corresponds to the Caroli-de Gennes-Matricon 
mode in the isotropic s-wave superconductors (in s-wave 
symmetry, J7o is 9 independent) . [j36|, |37j The quasipar- 
ticles with 9 away from the nodes in d-wave vortex are 
well localized near vortex cores and they are similar, in 
nature, with those in an s-wave vortex. As the angle 9 
approaches a nodal direction, however, the quasiparticles 
become more extended and farther away from the vortex 
cores. In this way the character of quasiparticles in d- 
wave vortex is very different from that of quasiparticles 
in the s-wave vortex. 

According to the theory by Kopnin and Volovik based 
on the relaxation time approximation, FFF resistivity is 
given by 



B 



Pf 



(n(9))T v n c 



(12) 



where (• • • ) denotes the average over the Fermi surface, 
r v is the relaxation time of quasiparticles, and n e is the 
carrier density in the vortex state. j26) In the theory of the 
relaxation time approximation, the transport coefficients 
are given in the form of the parallel circuit; the con- 
ductivity is expressed as a sum of the contribution from 
each part of the Fermi surface. Then magnitude of re- 
sistivity in vortex state expressed by Eq. ( |l2| ) is governed 
by the largest value of il(9) on the Fermi surface. This 
fact is physically interpreted in the following way. The 
quasiparticles with smaller fl(9) come from the vicinity 
of nodes. They are only weakly excited by vortex mo- 
tion, because such quasiparticles are extended in regions 
far away from vortex cores. On the other hand, quasi- 
particles with larger £1(9) are localized near vortex cores. 
Therefore it is likely that such quasiparticles are excited 
substantially by vortex motion and an appreciable devi- 
ation of the distribution function from the equilibrium 
state may occur. Thus when the gap has nodes, portions 
of the Fermi surface near the nodal directions do not 
contribute to (£1(6*)). This is in marked contrast to the 
isotropic s-wave superconductors, in which every part of 
the Fermi surface can contribute to (il(9)). The reduc- 
tion of the number of quasiparticles available for the en- 
ergy dissipation in the superconductors with nodes gives 
rise to the enhanced flux flow resistivity. This scenario 
has been adopted in Ref . Q to discuss the flux flow resis- 
tivity of UPt3. Although this argument explains the low 
field (H < 0.2H C 2) behavior expressed as (]?]), it gives 
no account for the v^ff-dependence of pj expressed as 
Eq.(^) observed in the almost whole regime at higher 
field (0.2id c2 < H < H c2 ). 

There is, however, another scenario. In the following 
part, we show that the reduction of r v in d-wave vor- 
tex states explains consistently both (Q) and (||). Here 
we regard the impurity scattering as the main process 
of relaxation in the cuprates. Within the Born approx- 
imation, r v is inversely proportional to the density of 
states (DOS) of quasiparticles available as the out-going 
states in the scattering process of localized quasiparti- 
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cles. On the other hand, the low energy DOS of quasi- 
particles in <i-wave vortex states are known to be larger 
than that in s-wave vortex states theoretically. [| , 26, 43] 
In Ref . |Bq| , Kopnin and Volovik calculated the density 
of states N V (E) per each d-wave vortex for energy E to 
obtain 



N V (E) ~ N o e (A/E) ~ N (r(E), 



(13) 



where Nq denotes the DOS on the Fermi surface in the 
normal states and r(E) = Hvf/E with vp the Fermi 
velocity. The singularity at E = is removed by a cut- 
off length. According to Ref. p(| , for energy E satis- 
fying r(E) > i?B with the intervortex distance Rb ~ 
^WH&jB, r(E) should be replaced by Rb for pure su- 
perconductors without impurity scattering. For impure 
and clean superconductors, instead, we speculate that 
r(E) should be replaced by Rb or the mean free path 
'v( = v f t v), whichever is smaller. We then expect that 



N v (0)/(N o e) 



ly/£: h < Rb 

\J B^jB, Rb < l v . 



(14) 



for E = 0. The quasiparticle DOS per each isotropic 
s-wave vortex is given by iv"o£o- fUl Therefore, the left- 
hand side in ([14]) gives the ratio of DOS in d-wave vortex 
states to that in the isotropic s-wave vortex state. From 
this fact and Eq. fl4|), we expect that 



t v (ci- wave) / t v (s-wave) 



£,/ly, lv < Rb 
\J B/H c2 , Rb < l v - 



(15) 



This reduction of the relaxation time in d-wave vortex 
also yields the enhancement of the flux flow resistivity 
Pf. If we assume here that this reduction of r v alone 
leads to the enhancement of pf, i. e. 

p/(d-wave)//9/(s-wave) ~ t v (s-wave) /r v (d- wave) (16) 

and pf (s-wave) ~ p n (B / H C 2), we obtain 

(l v /t)(B/H c2 ), l v <R B 



p f (d-wave)//9 n 



\/ B/H C 2, Rb < lv 



(17) 



We then see the upshot of the hyposesis ( [Tt| ) . The expres- 
sion ( fl7| ) is consistent with the experimental results on pf 
both in low fields Eq. (fjj) and in high field Eq. (||). From 



the relation l v ~ Rb at the crossover field 2T^3T from 
Eq. df) to Eq. (g), we obtain l v = 280 - 340A. From this 
value of l v and £ ~ 42 A (estimated from H C 2 = 20T), we 
obtain Z v /£ = 6.6 ~ 8. This value is somewhat larger 
than a ~ 2. With consideration of the crudeness of our 
estimation, however, we should say that these two values 
are of the same order. 

At the present state of the study, we do not know 
whether the dominant source for quasiparticle energy dis- 
sipation comes from the reduction of the number of the 
quasiparticles or the enhancement of the carrier scatter- 
ing rate. A detailed numerical calculation for the energy 
dissipation especially when each vortex overlaps with its 
neighborhood would be necessary. 

VI. SUMMARY 

The microwave surface impedance measurements in 
the vortex state of overdoped Bi:2201 demonstrate that 
the free flux flow resistivity in the moderately clean d- 
wave superconductor with gap nodes is remarkably dif- 
ferent from that in conventional fully gapped s-wave su- 
perconductors. At low fields, the free flux flow resis- 
tivity increase linearly with H with a coefficient which 
is far larger than that found in conventional s-wave su- 
perconductors. At higher fields, the flux flow resistivity 
increases in proportion to \J~H up to H C 2- Two possible 
scenarios are put forth for these field dependence; the 
enhancement of the quasiparticle relaxation rate and the 
reduction of the number of the quasiparticles participat- 
ing the energy dissipation in ci-wave vortex state. The 
present results indicates that the physical mechanism of 
the energy dissipation associated with the purely viscous 
motion of the vortices are sensitive to the symmetry of 
the pairing state. 
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